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METHOD AND APPARATUS FOR REDUCED 
RANK CHANNEL ESTIMATION IN A 
COMMUNICATIONS SYSTEM 



FIELD 

The present invention relates to wireless a 
More particularly, the present invention relates to a novel 
and improved method of reduced rank channel e; 
a communications system. 

BACKGROUND 



To improve the quality of wire 
munication systems often employ multiple radiating antenna 15 
elements at the transmitter to communicate information to a 
receiver. The receiver may then have one or more receiver 
antennas. Multiple antennas are desirable, as wireless com- 
munication systems tend to be interference-limited, and the 
use of multiple antenna elements reduces inter-symbol and 20 
co-channel interference introduced during modulation and 
transmission of radio signals, enhancing the quality of 
communications. The modeling, and thus design, of such a 
system, involves estimating several parameters of the space- 
time channel or link between the transmitter and receiver. 25 

The number of estimated channel parameters per 
transmit-receiver antenna pair is multiplied by the number of 
permutations of transmitter-receiver antenna pairs, creating 
increasingly complicated calculations and decreasing esti- 
mation quality. Therefore, it is desirable to have methods of 30 
channel estimation that use a reduced set of parameters. 
Similarly, there is a need for an improved method of channel 
estimation for radio communications systems having mul- 
tiple transmitter antennas. 

SUMMARY 

The presently disclosed embodiments are directed to a 
novel and improved method and apparatus for estimating 
channel parameters in a communication link in a wireless 40 
communication system having multiple transmitter antennas 
using a reduced rank estimation method. Each path from a 
transmitter antenna to the receiver constitutes a channel 
within the link. The number of channels, therefore, increases 
with the numbers of transmitter antennas and receiver anten- 45 
nas. The method exploits redundant and/or a priori knowl- 
edge within a system lo simplify the channel model used as 
a basis for the estimation calculations and to improve the 
estimation quality. In one embodiment, a covariance matrix 
is calculated and analyzed to determine if the number of 50 
channel parameters may be reduced for channel estimation. 
If not, all parameters are estimated, otherwise a reduced rank 
channel model is used for the calculation of channel param- 



represents a plurality of impulse responses between the 
transmitter and the receiver. Alternatively, the matrix may be 
a sample matrix describing the link. 

Further, determining the matrix may include estimating a 
plurality of parameters describing at least one channel. The 
parameters may include a distance between transmitter 
antennas. In one embodiment, the parameters include a 
transmittal angle with respect to a configuration ot the 
transmitter antennas. In an alternate embodiment, the matrix 
describes parametric relations of the link in the frequency 
domain. 

Further, ranking the matrix may include determining an 
eigenvalue for the matrix. In one embodiment, it the rank is 
equal to (NxM) a set of correlated impulse responses is 
applied for demodulating. In one aspect, a wireless appara- 
tus is operative to model a link in a wireless communication 
system by determining a matrix describing parametric rela- 
tions of the link; ranking the matrix; determining if the rank 
is less than NxM; if the rank is less then NxM performing 
an extraction of a subspace of the matrix; deriving channel 
impulse responses for each channel based on the extracted 
subspace of the matrix; and demodulating a received signal 
using the channel impulse responses. 

In another embodiment, a wireles 
ratus includes a correlator operative 1 
matrix representing a link with a transmitter based on signals 
received from the transmitter; a rank analysis unit coupled to 
the correlator and operative to estimate a rank of the 
covariance matrix; and a channel estimation unit coupled to 
the rank analysis unit and operative to generate a reduced 
rank channel estimate. The covariance matrix may represent 
a plurality of impulse responses between the apparatus and 
the transmitter. In one embodiment, the rank analysis unit is 
operative to determine an eigenvalue corresponding to the 
covariance matrix and is operative to compare the estimated 
rank of the covariance matrix to a predetermined full value. 

In still another embodiment, a method for ei 
link in a wireless communication system includes es 
a covariance matrix for the link; determining if the rank of 
the covariance matrix is reducible; reducing the rank of the 
covariance matrix; and estimating a set of impulse responses 
for the link using the reduced rank covariance matrix. 
Additionally, the method may include determining a corre- 
lation of the channel, ranking the covariance matrix; and 
performing an extraction of a reduced rank matrix out of the 



:t, a method for modeling a link in a wireless 55 
n system, the system having a transmitter 
having N antennas and a receiver having M antennas, each 
path from one of the N transmitter antennas to the M 
receiver antennas comprising a channel, includes determin- 
ing a matrix describing parametric relations of the link; 60 
ranking the matrix; determining it the rank is less than NxM; 
if the rank is less then NxM performing an extraction of a 
subspace of the matrix; deriving channel impulse responses 
for each channel based on the extracted subspace of the 
matrix; and demodulating a received signal using the chan- 65 
nel impulse responses. The matrix may be a covariance 
matrix describing the link, wherein the 



;less communication apparatus 
communication system having 
inas and a receiver having M 
; of the N transmitter antennas 
> comprising a channel. The 
of computer readable instruc- 
describing 



In one embodiment, a wir 
is operative within a wireless 
a transmitter having N ante 
antennas, each path from on 
to the M receiver antenna 
apparatus includes a first sel 
tions operative to determine 
the link; a second set of computer readable ii 
operative to rank the covariance matrix; a third set of 
computer readable instructions operative to determine if the 
rank is less than NxM; a fourth set of computer readable 
instructions operative to perform an extraction of a reduced 
rank matrix out of the covariance matrix if the rank is less 
then NxM; a fifth set of computer readable instructions 
operative to derive channel impulse responses for each 
channel based on the reduced rank covariance matrix; a sixth 
set of computer readable instructions operative to demodu- 
late a received signal using the channel impulse responses. 
The apparatus may further include an equalizer operative in 
response to the sixth set of computer readable ir 
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wherein a configuration of the equalizer is determined by the 
rank of the covariance matrix. In one embodiment, the 
apparatus includes a seventh set of computer readable 
instructions operative to derive a correlated channel impulse 

In still another aspect, a wireless communication appara- 
tus includes a channel estimation means operative to esti- 
mate a covariance matrix representing a link with a trans- 
mitter based on signals received from the transmitter; a rank 
analysis unit coupled to the correlator and operative to 1Q 
estimate the rank of the covariance matrix; and a channel 
estimation means coupled to the rank analysis unit and 
operative to generate a reduced rank channel estimate. 

Further in another aspect, a wireless communication 
apparatus includes a correlator operative to estimate a cova- ^ 
riance matrix representing a link with a transmitter based on 
signals received from the transmitter; a rank analysis unit 
coupled to the correlator and operative to estimate the rank 
of the covariance matrix; and a channel estimation means 
coupled to the rank analysis unit and operative to generate . 
a reduced rank channel estimate. 

In yet another aspect, a method for estimating a link in a 
wireless communication system includes estimating a cova- 
riance matrix for the link; determining if the rank of the 
covariance matrix is reducible; reducing the rank of the ,, 5 
covariance matrix; and estimating a set of impulse responses 
for the link using the reduced rank covariance matrix. The 
method may further include determining a correlation of the 
channel, ranking the covariance matrix; and performing an 
extraction of a reduced rank matrix out of the covariance , Q 

In another embodiment, a wireless apparatus includes 
channel estimation means operative to determine significant 
delays and determine a set of estimates of full dimension 
channel parameters associated with the significant delays, 35 
wherein each one of the set of estimates corresponds to an 
instance in time; eigenvalue computation means operative to 
determine eigenvalues of the set of estimates of the full 
dimension channel parameters and find any dominant eigen- 
values; and channel estimation means operative to determine 40 
a set of reduced rank channel parameter estimates in 
response to the dominant eigenvalues. Further, the apparatus 
may include eigenvector computation means operative to 
determine at least one eigenvector associated with one of the 
dominant eigenvalues of the set of estimates; wherein the 45 
channel estimation means uses the at least one eigenvector 
to project the set of estimates of the full dimension channel 
parameters onto the subspace spanned by the at least one 
eigenvector. 

BRIEF DESCRIPTION OF THE DRAWINGS 50 

The features, objects, and advantages of the present 
invention will become more apparent from the detailed 
description set forth below when taken in conjunction with 
the drawings in which like reference characters identify 
correspondingly throughout and wherein: 

FIG. 1 illustrates configurations of wireless communica- 
tion systems including multiple transmitter antennas; 

FIG. 2 illustrates a model of a wireless communication 
system according to one embodiment; 6Q 

FIG, 3 illustrates a model of a channel between transmit- 
ter and receiver in a wireless communication system; 

FIG. 4 illustrates the physical layout of antennas in a 
transmitter of a wireless communication system; 

FIG. 5 illustrates a flow diagram of a method of reduced 65 
rank channel estimation for a wireless communication sys- 
tem according to one embodiment; 
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FIG. 6 illustrates a plot of the estimation gain of one 
embodiment; and 

FIG. 7 illustrates a system configuration according to one 
embodiment. 

FIG. 8 illustrates an exemplary embodiment of a wireless 
communication system. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

Multiple radiating antennas may be used to improve 
transmission quality in a wireless communications system. 
In the design of third generation mobile radio systems, for 
example, various transmitter antenna diversity techniques 
are presented. Multiple transmitter antennas may be used to 
communicate information to a receiver using a single or 
multiple receiver antenna(s). Multiple antenna systems offer 
an improvement in quality. However, the improvement is 
dependent on the accuracy of the channel model used in the 
receiver to demodulate the transmitted information. Model- 
ing oi the transmission channel uses parameter estimates and 
determines an effective channel impulse response tor the 
channel. When multiple antennas are used, the modeling 
involves estimates of each transmission channel for all 

The transmission channel from transmitter to receiver is a 
space-time channel described generally by at least one 
impulse response. Often there is little change in the channel 
parameters from one channel to another, such as where the 
channel impulse responses differ only in phase. In such a 
case, it may not be necessary to derive estimates of impulse 
responses independently for each channel, but rather some 
information may be reused. When channels are correlated, a 
reduced rank representation of the channels may be used. 
Reduced rank refers to the reduced number of completely 
uncorrelated channels used to describe the link between 
transmitter and receiver. One way to observe this reduced 
rank is the rank reduction of the channel covariance matrix 
used to describe the mutual statistical dependencies of the 
different channel impulse responses. Note that the reduced 
rank can also be realized by other parameter measures. For 
example, in one embodiment a sample matrix is formed of 
columns comprising samples of channel impulse response 
estimates over time, wherein the reduced row rank of such 
a sample matrix is applied as described herein. A reduction 
in rank may result in a less complex filter or demodulator, 
i.e., reduces the number of filters and/or filter elements 
and/or demodulation units used in the receiver. Furthermore, 
the reduction of the number of estimated parameters used to 
characterize the channel leads to improved accuracy of the 
channel model. 

FIG. 1 illustrates configurations for wireless communica- 
tion systems having multiple transmitter Tx antennas. Two 
paths are illustrated: a first multiple input, multiple output 
(MIMO) and a second path multiple input single output 
(MISO). The MISO configuration places multiple Tx anten- 
nas in communication with a single Rx antenna. The MIMO 
configuration extends this to multiple Rx antennas. The 
channel model for one of the systems of FIG. 1 is illustrated 
in FIG. 2, in accordance with one embodiment, specifically, 
for a wireless system employing coherent demodulation, 
having a link between a transmitter and a receiver, wherein 
at least the transmitter employs multiple antennas. The 
wireless communication system 10 includes a transmitter 12 
and receiver 16 that communicate via an air interface. A 
channel model 14 represents the channels for antenna pairs 
between transmitter 12 and receiver 16. Channel model 14 
considers the channels within a link, such as the MISO link 
of FIG. 1. 
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Continuing with FIG. 2, let N Tx be the number of antennas 
used at the transmitter 12 and N^. the number at the receiver 
16, respectively. In general, for each significant propagation 
delay between "transmitter and receiver, (N n -N & ) transmis- 
sion channels exist for the pair, wherein for a significant 5 
propagation delay the received signals resemble the known 
transmitted signals with high certainty. In other words, 
define N £ as the number of significant propagation delays, 
also referred to as echoes. The (N^-N^-N^) channel 
impulse response samples are then estimated to perform JQ 
coherent demodulation. When the channels are uncorrelated, 
the (N^/Njjj.-N^) channel impulse response samples are 
modeled as completely uncorrelated random processes and 
the estimates of these channel impulse response samples 
may be derived independently without loss of demodulation 
performance, However, if the fN^-N^-NJ channel impulse 15 
response samples are not uncorrelated random processes the 
(N-^-N^'N^) channel impulse response samples may be 
modeled as a linear combination of a smaller number N c ,, of 
channel impulse response samples, wherein N c/l <(N :r _ > .'N Kl . 
N E \. Such cases include, but are not limited to, minimal 20 
angular spread at the transmitter and/or receiver in the 
effective channels due to propagation conditions. If N CA is 
known, or estimated, and the linear transformation of the 
N cfl channel impulse response samples are resolved into the 
(N^-N^'N^) channel impulse response samples, then mod- 2 s 
eling may be accomplished with the N ch channel impulse 
response sample estimations. This reduces the number of 
parameters to be estimated while increasing estimation 
quality, yielding an increase in the demodulation perfor- 
mance. Even if the exact representation of the linear trans- 30 
formation of the N cfe channel impulse response samples into 
the corresponding (Nj^-N^-N^) channel impulse response 
samples is not known, the modeling may still be accom- 
plished with channel impulse response sample estima- 
tions if the subspace spanned by vectors of this linear 35 
transformation is known or can be estimated. 

This principle is referred to as "Reduced Rank Channel 
Estimation." The transformation of the N c/l uncorrelated 
channel impulse responses into Ihe (N^-N^-N^) correlated 
channel impulse responses can depend on factors including. 40 
but not limited to, the antenna configuration, antenna 
patterns, polarization characteristics, propagation conditions 
and more. In some cases the transformation might be known 
a priori, in other cases it can be derived or estimated, for 
example by anglc-of-arrival estimation. The subspace 45 
spanned by the linear transformation of the N ch channel 
impulse response samples into the corresponding 
(N ZT N RT N f ) channel impulse response samples can be 
determined by estimating the rank and the eigenvectors of 
the (N^-N^'N^-dimensional covanance matrix of the 50 
(N^'N^'N^) channel impulse response samples. This sub- 
space can also be determined by using a singular value 
decomposition of a matrix holding columns with all 
(N & -N Rt N £ ) channel impulse response sample estimates 
for different points in time. Note that if the channel impulse ss 
response samples are corrupted by a known correlation 
noise, and if the noise correlation can be estimated, the 
(N^-N^lSy channel impulse response samples may be 
filtered by a noise de-correlation filter first. 

In one embodiment, the rank reducing transformation is 60 
known a priori or is estimated. In other words, the mapping 
of the N cjl channels onto the (N^-N^) channels is ascer- 
tainable. The reduced rank channel is then estimated using 
the ascertained transformation. When desired, an equivalent 
full dimensional channel model may then be derived from 65 
the reduced rank estimate by transforming the reduced rank 
estimate back to the larger dimension. 



In an alternate embodiment, the rank reducing transfor- 
mation is not directly known, but the subspace spanned by 
the transformation may be extracted from the dominant 
eigenvectors of the channel covariance matrix. Note that the 
subspace may be referred to as the signal subspace or the 
channel subspace. The process involves first estimating a 
channel covariance matrix and finding the dominant eigen- 
values. By determining the associated eigenvectors which 
span the channel subspace, the process projects the conven- 
tional channel estimate into the channel subspace, yielding 
a reduced rank channel model with reduced estimation 
errors. If desired, the reduced rank model may be trans- 
formed back into an equivalent full dimension channel 

FIG. 3 illustrates a model 18 of a MIMO channel for 
continuous time having a linear MIMO filter 20 with N & 
inputs and N Kt outputs. The linear MIMO filter 20 is defined 
by the (N-^xN^,) matrix H(t) comprising of linear functions 
h,/t), i=l . . . N Tj , j=l . . . IV Generally, h- (t), i-1 . . . N^, 
j=l . . . Njfr are unknown linear functions. The linear MIMO 
filter 20 represents the (N rx -N^) radio channels through 
which the N rT transmit signals pass to the N Sr receiver 
antennas. These radio channels are characterized by their 
channel impulse responses h,-,{t), i-1 . . . N & , j=l . . . fV 
The input signal to the model, x (t), is a (N r:t xl) column 
vector representing the N rx band-limited transmit signals, 
and the output signal from the model, y (t), is a (N fc xl) 
column vector, sampled at t=T, 2T . . . , as illustrated by- 
switch T, where the bandwidth of the transmitted signals is 
less or equal to 1,'T. The received signals contain additive 
perturbation signals represented by the (N^xl) column 
vector l (t), introduced due to noise or co-channel interfer- 
ence. The additive perturbation signals are added at sum- 
mation nodes 22. The relation between the input signals 
x ( 0, the channels H(t), the perturbation z (t) and the output 
signals y (t) is given by 



y (i)=H T (t) 



(1) 



where * denotes the convolution. 

FIG. 4 illustrates the physical configuration of antennas at 
the transmitter of an exemplary embodiment modeled as in 
FIG. 2. A reduced rank method is applied to estimate the link 
represented by channel model 14, having a transmitter 12 
configured with the four (4) antennas, each spaced at a 
distance "d." The specifics of the configuration and model 
are discussed hereinbelow. Note that the estimation proce- 
dure is performed at the receiver 16. A reference direction is 
given by the horizontal line. Angles of transmission are 
measured with respect to this reference. The angle "a" 
corresponds to an angle of a propagation path with respect 
to the reference within a 2-D plane as illustrated. A range of 
angles with respect to the reference is also illustrated. The 
following method is used at the transmitter 12 in system 10 
to estimate the link. 

FIG. 5 illustrates a flow diagram of an exemplary method 
of channel estimation used to process signals in a receiver 
unit in accordance with one embodiment. Process flow 
begins by searching for significant propagation delays in the 
channel, i.e. searching for significant echoes at step 40. In 
one embodiment the process involves a sliding correlation of 
the received signals with known transmitted signals or 
known components of the transmitted signals. Correlation 
refers to the degree with which the received signals are 
related to the known transmitted signals, wherein a perfect 
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correlation proves a relationship between the signals with 
high confidence. For time-shitted signals, wherein sliding 
delays are used to shift the received signals in time, a 
resultant sliding correlation provides the degree of certainty 
with which the time-shifted signals resemble the known ; 
transmitted signals. Thus in the wireless system context, 
sliding correlation relates to the synchronization of known 
signals transmitted by the Tx antennas with time-shifted 
versions of the received signals. The exemplary embodiment 
of reduced rank channel estimation uses sliding correlation 
of the received signals with known transmitted signals to 
estimate the number N £ and the values x lt x 2 . . . ,x N of 
significant propagation delays, i.e. delays for which the 
received signals shifted back by these delays in time ^ 
resemble the known transmitted signals with high certainty. 
The procedure of sliding correlation in order to find signifi- 
cant propagation delays is also known as "searching" in 
CDMA systems. 

The method then estimates parameters for multiple 2' 
observable channels between the N rv transmitter antennas 
and the receiver antennas at step 42. The channels are 
radio network connection pairs coupling at least a portion of 
the N n . transmitter antennas to at least a portion of the 
receiver antennas. In the exemplary embodiment, there is a 2. 
connection between each transmitter 12 antenna and each 
receiver IS antenna, resulting in (N^-N^) channels. The 
parameters describing the multiple channels are those char- 
acteristics that impact the impulse responses of the channels. 
Assuming that N £ significant propagation delays (echoes) 3' 
exist between transmitter and receiver, (N^-N^N^) com- 
plex samples of the (Nj^-N^) channel impulse responses 
could be used as a set of parameters describing the multiple 
channels. This set of parameters is denoted by a 
((N^N^-N^xl) vector termed ~h herein. The relation 3 
between x (t), h , z (t), and y (t) is developed hereinbelow. 

With x±, x 2 , . . . ,x N being the significant propagation 
delays between transmitter and receiver, the model 
described by (f ) can be expressed as 4> 



v(0 = 5>V.)-*r-, 



Since the output signals y (t) are sampled at a sampling rate 
of 1/T, vectors containing the discrete time samples can 
represent segments of a finite duration of the continuous 
time signais. For the sake of simplicity, the received signals 
y (t) are described herein by a discrete time representation 
over a finite duration of time t=0, T, . . . (N^l)T, where N r 
is the number of samples taken over time. Therefore, the 
toflowing abbreviations are used. Each discrete time trans- 
mitted signal at antenna n delayed by x is given by a vector 



'>M4-v:c- 



xJT-V 



-^((Nj-rir-Ttr, , 



Wherein the matrix describing all of the discrete time 
transmitted signals delayed by x is given as 



The matrix A describes all discrete ti 
having significant delays. 



ar describing each discrete t: 



..(/<«.)© S(T„ E ))] (9) 

ic perturbation signal at 



„=[„< 



(II) 



The vf 



ir of the disc 



lived signal at ai 



This can be transformed ir 



is given as 

T^MyMrJP) ■ yJ&r-WT, ■ ■ ■ (12) 
and the vector of all discrete time received signals is given 



yCO-TP*' ©T^-tj) (/(*»© T r ( 
? r ('-% £ ))]-7+?(0, 
where (x) denotes the Kronecker tenso 
(NjkXNjjJ identity matrix, and the ve 
from the matrix H(t) such that [[~h ]] 



Using the above abbreviations, the discrete time output 
signals of the MIMO channel model 18 illustrated in FIG. 3 
over a period of time from t=0,T . . . , (N r -f)T may be 
reduced to the simple model 



S=[h,„K) 



hold. The ((N^-N^xt) vector h , otal (x e ) is consisting of the 
elements of the matrix H(t) sampled at x c with all columns 



of H(T e ) stacked on top of each other ir 
which is denoted by the operator ve 



the vector h rot Jx e ), 
i{U(x e )} in (5), i.e. 



The second step in the flow diagram in FIG. 5 at step 42, is 
to repeatedly process estimates for a set of parameters 
characterizing the multiple channels between transmitter 
and receiver. For the above-described mathematical repre- 
sentation of the channel model, this may be equivalent to 

processing estimates h t " ) , n=f . . . N,, of the vector h in 
(14) for different points in time. A conventional method 
> uses the correlation of the received signals, shifted back in 
time by certain delays, with known transmitted signals, such 
as pilot signals specific to 
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predetermined training sequences. As the significant propa- 
gation delays t 1 ,t 2 , . . . ,x N are already determined in step 
40, the exemplary embodiment of reduced rank channel 
estimation uses the correlation of known transmitted signals 
with versions of the received signals, shifted back in time by ; 
t-L.tj, . . . ,t Vt , to generate a channel model, such as channel 
model 14 of FIG. 2, characterized by the vector h . If the 
noise vector n represents spatial and temporal white 
perturbation, wherein the noise covariance matrix is given as \ , 
R„=< n n H >=a l -V f '^'' vr >, and if the matrix A comprises of 
the a priori known signals, such as pilot symbols of a CDMA 
system, channel estimates obtained by correlation can be 
described by 



If the noise vector n does not represent spatial and temporal 
white perturbation, the channel estimates obtained by cor- 
relation can be described by : 



Note that R„ might be known a priori or could be estimated 
from the received signals. The channel estimate of (f6) 2 
contains a perturbation vector A^R," 1 n with a covariance 
matrix of R^=A H R„- 1 A. This covariance matrix is not 
diagonal in general, i.e., the components of the perturbation 
vector contained in h are correlated in general. If R p is 3 
known or can be estimated, the components of the pertur- 
bation vector contained in h could be de-correlated by 
transforming h with R^f 1 ' 2 . This will be assumed in what 
follows, wherein 3 



shall hold. 



g is a (N CA xl) vector with uncorrelated components, i.i 
R g =<"g • "g"> is a diagonal (N CA xN cfc ) matrix. In this casi 
the channel covariance matrix is given as 



R h =BRB". (20) 

As a consequence, the rank of the channel covariance 
matrix R h is equal to N c/r Given (20), and assuming that 
correlation according to (17) is used to derive the channel 
impulse response estimates h , the covariance matrix of 
h is given as 

s R k =<t-t^R n v2h B-R e B H R/ !2 +f XT '-- Wli ' Nl: \ (21) 

Due to the reduced rank N ch of R h , the eigenvalue 
decomposition 

R p im B-R s -B H R p v2 =R p v2,i R h R jl v2 =EKE H , (22) 

yields only N ch non-zero eigenvalues, where A is a diagonal 
matrix containing the eigenvalues and E is a square matrix 
containing the eigenvectors of R p 1/2ff B-R g B' 7 R £ 1 '' 2 . With 
(21) and (22) the covariance matrix estimate R,, maybe 
5 expressed by 



i.e., R,, shares the eigenvectors with R p 1!2H -K-R g -B H R p ln . 
Since A is a diagonal matrix with only N cv , non-zero 
elements, (N T:) .-N R -N E )-N C ^ eigenvalues of R h are constant, 
and N c/l eigenvalues of R A are larger than the former ones. 
These larger eigenvalues are termed dominant eigenvalues 
in the sequel. With a diagonal matrix A c , containing all 
dominant eigenvalues of the estimated channel covariance 
matrix, the matrix E c , containing the corresponding 
eigenvectors, and with the matrix E N , containing the remain- 
ing eigenvectors, (23) becomes 



R h -E c 



k of the 



As illustrated in 
channel parameters is estimated at step 44. Covariance 
measures the variance of one random variable with respect 
to another. In this case, the covariance matrix describes the 
variance of the various channel parameters with respect to 
each other. According to the above -described mathematical 
representation of the channel model, step 44 corresponds to 
processing an estimate R h of the channel covariance matrix 
R 4 =<~h ■1Th>. Such 



e may be given a; 



If the MIMO channel has a reduced rank wherein N c/l < 
(N^-N^'N^), i.e. the (N^-N^-NJ MIMO channel impulse 
response samples can be described as a linear combination 
of N CA uncorrelated channel impulse response samples. The 
channel vector h can be modeled as a linear transformation 



of a channel vector g of reduced dim 



Therefore, the matrix E c contains the eigenvectors spanning 
the channel or signal subspace. 

The estimated covariance matrix R,, is then ranked at step 
46, meaning that the number of dominant eigenvalues is 

45 estimated. The rank is compared to a maximum value 
"MAX" at step 48, MAX is equal to the total number of 
estimated channel parameters in the vector h . In other 
words, MAX is equal to (N^ N^-N^). As many of the 
mechanisms impacting correlation, such as the directionality 

50 of the propagation paths, do not change quickly over time, 
the correlation characteristics may be estimated by averag- 
ing over rather long time intervals in comparison to the 
inverse fading rate of the channel(s). 
The rank of the covariance matrix determines whether the 

55 (Nji'Nfc'r^) channel parameters describing the (N^-N^) 
existing transmission channels can be modeled as a linear 
combination of a smaller number N ch of equivalent uncor- 
related channel parameters. If a reduced rank is available, 
the channel subspace E c of the estimated covariance matrix 

60 R h is derived at step 52. Note that instead of using the 
estimated covariance matrix R h , the rank of R h and the 
channel subspace E c can also be derived from the matrix of 
channel parameter e: 



and wherein B is a ((N^-N^'N^xN^) matrix describing 
the linear transformation. As given hereinabove, the vector 



(25) 



by using singular value decompositioi 
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With the chanr 
nel parameter vet 



1 subspace E c , reduced dimension chan 
>rs are estimated at step 54, according tc 



effectively projecting the originally estimated channel 
parameters into the channel subspace. This projection into 
the channel subspace reduces the estimation error. If a 
reduced complexity demodulator is used in the receiver, 
which uses the reduced rank channel, i.e., takes only a 
reduced number of channel parameters into account for 
demodulation, the estimates of (26) may be directly used in 
the demodulator for coherent demodulation. In other words, 
processing would flow directly from step 54 to step 58, or at 
a minimum step 58 would used the reduced rank estimates. 

If a conventional receiver, designed for the full rank 
channel model, is to be used, the estimates g ^ may be 
transformed back into the full dimensional space at step 56, 
according to 



wherein the factor R p " 1/2/f is used to make the estimate 
unbiased. Note that the estimate of the channel subspace E c 
may be updated continuously by using a sliding time win- 2 
dow for the estimates R,, or X ; , respectively. This eliminates 
the delay of waiting for a new complete sample set, by using 
a portion of the previous sample set with incrementally 
time-shifted new values. 

If rank reduction is not possible, processing continues to 3 
use the full rank of the system to model the channel at step 
50. In this case the method estimates the (N rx N KT N £ ) 
channel parameters independently from each other. Once the 
system is modeled, signal demodulation continues at step 
58. 3 

The MISO path illustrated in FIG. 1 is provided as an 
exemplary embodiment. As illustrated, the transmitter, Tx, 
has four (4) radiating antennas (N & =4) and the receiver, Rx, 
has one (f) antenna (N^-i). For simplicity, several assump- 
tions allow a straightforward analysis demonstrating the 4 
applicability of the exemplary embodiment to modeling a 
system as illustrated in FIG. 1. First, the example assumes 
that each 'Tx antenna transmits a pilot signal specific to that 
antenna, wherein the antenna-specific pilot signal is time- 
aligned and orthogonal to the pilot signals of the other Tx 4 

Second, assume the channels are frequency non-selective 
fading channels, each made up of a large number, P, of radio 
network paths. The paths each have approximately a same 
run length and a same attenuation. The second assumption 5 
ensures that the relative propagation delay is smaller than 
the inverse of the transmission bandwidth. The propagation 
delay of two radio paths is typically due to differences in run 

Third, the channel model is restricted to 2-D propagation, 5 
i.e. all effective radio paths are located in a 2-D plane. See 
FIG. 4. Additionally, the geometry of the effective radio 

wherein each path departure angle, measured with respect to 
a reference direction of Tx, are concentrated around an 6 
average angle, a. The radio path angles are Gaussian dis- 
tributed having mean a and a standard deviation a. For one 
simulation, a is selected randomly between -60 and +60 
degrees. The standard deviation a is assumed to be square 
root of two degrees. Fourth, the arrival paths at Rx are 6 
assumed uniformly distributed between 0 and 360 degrees to 
consider local scattering. Fifth, no line of sight exists. 



Sixth, assume a specific phase and Doppler shift for each 
path. The path-specific phase is selected randomly according 
to a uniform distribution between 0 and 2it. Additionally, the 
path-specific phase is adjusted for each Tx antenna accord- 

5 ing to the geometrical antenna configuration, i.e., the 
antenna location with respect to a reference point. For phase 
adjustment, assume object scattering is considered in the far 
field. The channel-specific Doppler shifts are generated 
according to a uniform distribution of the angles of arrival 
paths at Rx, a carrier frequency and a predetermined Rx 
speed. In the exemplary embodiment, the carrier frequency- 
is assumed to be t.8 GHz and the receiver speed equal to 60 
km/li yielding a maximum Doppler shift of f 00 Hz. In the 

jj exempfary embodiment, each Tx antenna covers a 120- 
degrec sector, with the antenna patterns all oriented towards 
a=0. 

Given the exemplary system as detailed, application of 
the process of FIG. 5 provides a channel model having a 
20 time variance according to the classic Doppler spectrum. It 
is possible to consider an antenna-specific radiation pattern. 
With this channel model, the channel impulse responses for 
the channels seen through the different transmitter antennas 
can be generated using the same set of radio paths, thus, 
introducing realistic correlation in the fading of the different 
channels. 

On the receiving side of the air-interface, at the single 
antenna of Rx, the method derives an impulse response 
3 0 estimate for each of the four transmission channels, i.e., the 
four radio network connections between Tx antennas and the 
Rx antenna. The estimate is based on the a-priori knowledge 
of spreading codes used to generate the antenna-specific 
pilot signals associated with each Tx antenna. 
35 Referring again to FIG. 4, in the geographical configura- 
tion of antennas at Tx., the antennas are positioned in a line 
having constant spacing d between neighboring antennas, 
wherein d=X, i.e., antennas are spaced one wavelength apart. 
Note that Rx has a single omnidirectional antenna. A total 
40 number of effective radio paths is considered with P=50. 
Channel specific variabies, a , f , and $ , represent, 
respectively, the angle measured from the reference line, the 
Doppler shift and the phase. The equation describing the 
channel impulse response for a channel between Tx antenna 
45 n and the Rx antenna is given as 



("- K 



wherein g„(a) is the antenna-specific complex azimuth 
radiation pattern of each Tx antenna. 

If a channel has no angular spread, and all path-specific 
angles a p are equal to a, the channel impulse response for 
each Tx antenna is given as 
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-continued 



which the pilot signals a 
following holds 



orthogonal, and wherein the 



where h(t) is the equivalent channel impulse response for an 
equivalent isotropic Tx antenna at the reference point. 

In this case the channel impulse responses for the different 
Tx antennas only differ by a complex factor, i.e., the 1 
channels are completely correlated. The steering vector is 
then defined as 



and the receiver noise signal z(t) represents white Gaussian 
noise. The signal received by the single Rx antenna is 
described as 



and the channel impulse response v 



The four (4) channel impulse responses seen from the Tx 
antennas are then copies of the channel impulse response 
h(t), weighted by four different complex factors, which 
means, the vector h (t) is a linear transformation of the 3 
scalar h(l) given by 

i.e., the vector g in the linear transformation of (ly) is in 
this example equal to the scalar h(t) and the matrix B is equal 
to the vector t*(a). This means the channel covariance 
matrix R A -<7-"?*) is equal to R ; =T(a)T(af <|h(t)| 2 > in 4 
this example. If the steering vector a (a) is known, such as 
a-priori knowledge of the antenna configuration and the 
radio path direction a, it is sufficient to estimate the scalar 
h(t) and either calculate an estimate for h (t) using the linear s 
transformation with a (a) or use the estimate of h(t) and a( 
a) directly for demodulation. 

Note that for the case when a (a) is known it may be 
sufficient to estimate h(t) and then compute an estimate of , 
h (t) from the scalar estimate of h(t). If the demodulator is 
designed such that the channel consists of a single scalar, 
i.e., the demodulation considers a (a), then it is possible to 
demodulate using a (a) and the scalar channel. 6 

The antenna-specific pilot signals at the transmitter are 
termed x„(t), n=l . . . N Tx , and the relationship is defined by 



Conventionally, correlating the received signal with the 
four (4) pilot sequences derives a set of four (4) channel 
estimates. Wherein the pilot signals are orthogonal over a 
pilot symbol period, this estimation is then repeated at the 
pilot symbol rate. Such a correlation procedure is generally 
, referred to as "integrate and dump" and may be expressed as 

(37) 



wherein h ca J"> is a vector made up of cc 
integrate and dump, channel estimates derived from the n-th 
pilot symbol. If (34) is transformed into a discrete time 
representation, by putting N r =T/T samples of the pilot 
signals x„(t) into the columns of the matrix A, NT samples 
of the noise signal z(t) into the vector n , and N r samples of 



the received signal y(t) into the v 



r , (34) yields 



..N Tx y 



Then the discrete time representation of (37) is 



If the channel v 
neglected, (39) becomes 



Considering the linear transformation of h(t) in 
received signal is expressed as 



The pilot signals are made up of segments, each having a 
duration J s , referred to as the pilot symbol duration, over 
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From this, an estimate of the scalar h(t) is derived as 




Again, when the channel variations within one pilot symbol 
are ignored, (13) becomes 

* " i ' ( " ri)+ ,v r .||H(s)|| ! A 

From this scalar estimate, using the linear transformation, a 
new estimate of the channel impulse vector is generated as 




Ignoring channel variations within one pilot symbol, (44) 
becomes 




(45) 



If a (a) is not known a-pnori, it may be estimated using 
the covariance matrix given by 

R„=<TlnT s )-h' T {nT s !>=7(a)-7' T (aYP,,, (46) 

with V h being the average power ot the scalar channel 
impulse response h(t). The covariance matrix R A can be 
approximated as 

V, (47) 



7,270 Bl 

16 




Note that the vector v max in this example is equal to the 
channel subspace E c . In general, the estimate R h is used to 
determine whether the rank of the channel estimation cova- 
riance matrix can be reduced. If R h is full rank, the channel 

10 estimation problem is not reduced to a smaller dimension. 
According to the exemplary embodiment, orthogonal 
pilot signals of binary chips have a chip rate of 1 .2288 Mcps, 
and a pilot symbol duration of 64 chips. With this channel 
model, a received signal, including white Gaussian noise, is 

15 generated for 4000 consecutive pilot symbols having a pilot 
Signal-to-Noise Ratio (SNR). From the received signal, 
4000 conventional vector estimates, h are generated. 

The thus generated covariance matrix is averaged over 
20 these 4000 consecutive conventional channel estimates. In 
the exemplary embodiment, the process takes approximately 
208.3 ms. After extracting the eigenvector corresponding to 
the maximal eigenvalue of R A , the matrix R 0 is calculated. 
Subsequently, 4000 new vector estimates h mr J-"' > are pro- 
duced according to 

"VT„„/")=jj J -t £OTF ; " ) . (50) 
Using the exemplary embodiment, iterations are repeated 
N , 7 =50 times. Ovcr_t_hc 50 iterations the transmitter angles 
are varied such that a is uniformly distributed within (+/- 
60) degrees, while the angular spread remains constant, 
having a standard deviation of square root of two (V2) 
degrees. Additionally, the channel parameters for a given 
pilot SNR are varied. The varied parameters represent radio 
path direction(s), path-specific phase, and path-specific 
Doppler shift, for a certain pilot SNR. An equal number of 
iterations is performed for different pilot SNR values. A 
comparison of the quality of the set of conventional esti- 
4Q mates to the set of new vector estimates, with respect to the 
reduction factor of the mean squared estimation error that is 
averaged over time and iterations, is made using the esti- 
mation gain given as 



which averages the vector with the conventional channel 
impulse response estimates over a number, N Jvm , of pilot 
symbols. 

For the case without noise and having an angular spread 
equal to zero, R h is rank one (f) and the vector a (a) spans 
R h . Thus (47) reduces to 




Note that the normalized vector a (ci)/|| a (a)|| spans R a . 

For a noisy case with sufficient low noise power and 
sufficient low angular spread, R h is still dominated by one 
eigenvalue. Therefore, the process performs an eigenvalue 
decomposition of R ; ,. When one eigenvalue is much larger 
than all other eigenvalues, it is an indication that the angular 
spread around a (a) was rather small. Therefore, as v ,„ OT is 
the eigenvector corresponding to the largest eigenvalue of 
R h , the approximation becomes 



45 




For the exemplary embodiment, FIG. 6 illustrates the esti- 
mation gain in dB as a function of the pilot SNR. Wherein 
the pilot SNR is defined as the ratio of the average energy 
per pilot chip E c of one pilot signal received at the single 
^ antenna receiver to the received noise power density I Q in 
dB. 

The upper limit for the estimation gain is determined by 
the number of transmit antennas, which is illustrated in FIG. 
6 as 6 dB. As illustrated in FIG. 6, the estimation gain 
approaches the upper limit even though the assumed angular 
spread is not zero and the received signal is severely 
corrupted by noise. The reduction of the estimation gain 
with increasing pilot SNR is due to the non-zero angular 
spread. 

Although the channel impulse responses are not com- A 
pletely correlated, the derivation of the impulse response 
h assumes this property. For larger angular spreads, a 
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smaller estimation gain is expected. For small angular 
spreads, the estimation gain appears considerable. Note that 
in general, for residential and suburban environments, a 
standard deviation of one (1) to two (2) degrees is frequently 
observed. Note also that it is possible to evaluate the 5 
performance improvement of the reduced rank channel 
estimation method using a Monte-Carlo-simulation to derive 
the reduction of estimation errors of the channel impulse 
responses as compared to conventional channel estimation 
using independent correlators. to 

Reduced rank channel estimation for systems using mul- 
tiple transmitter antennas allows improvement of the chan- 
nel estimation quality under certain propagation conditions 
with limited diversity due to correlated fading. As the 
mechanisms affecting correlation, such as the directionality 15 
of the radio wave propagation, change relatively slowly over 
time, the correlation characteristics may be estimated by 
averaging over extended time intervals. This is in contrast to 
the time intervals associated with inverse fading rate of the 
channel and thus allows improved accuracy in estimating the 20 
correlation characteristics. 

Reduced rank channel estimation for multiple transmitter 
antennas is also applicable to frequency-selective channels 
by computing either separate estimates of the correlation 
characteristics or by computing estimates of the correlation 25 
characteristics across all propagation delays. Separate esti- 
mates refers to computation of R A , for each propagation 
delay. Reduced rank channel estimation is then performed 
taking rnto account each delay occurring in the trequency- 
selecuve channel impulse response. In an alternate 30 
embodiment, wherein additional information, such as the 
antenna conhsuratron at the transmitter, is known a-pnorr. 
the step ot estimating the linear transformation ot the 
reduced number ot uncorrected channels into the larger 
number ot correlated channels may be more accurate. 35 
Additionally, the reduced rank estimation process may be 
extended to cases with more than one rcccrvcr antenna. In 
this ease, the estimation is performed lor the MIMO 
channels, as illustrated in Fib. 1. While the present example 
involves a system employing coherent demodulation. 40 
reduced rank channel estimation as descrrbed herein is also 
applicable to communication systems employing non- 
coherent demodulatron. 

A receiver 100 according to one embodiment ot the 
present invention is illustrated in FIG. 7. The receiver 100 45 
has a simile antenna 102 that receives signals irom a 
transmitter having multiple antennas. ITie received sisnals 
are first processed by the preprocessor 104. The signals are 
then provrded to a correlator 106. which is used as a shdinn 
correlator for searching and as a correlator for the significant 50 
delays for channel estimation. In an alternate embodiment 
the delays are determined in software without use of a 
correlator. The outputs of the correlator 106 are used to 
provide an estimate of the covariance matrix. In one 
embodiment, the correlator 106 is made up of fingers to form 55 
a rake, having one finger for each combination of transmitter 
antenna, receiver antenna and significant delay. The esti- 
mates are provided to the central processor 112 via bus 116. 
The processor 112 stores the channel parameter estimates in 
memory 114 so that the estimates may be used to derive the 60 
channel covariance matrix averaged over time. 

From memory 114, the estimated 
provided to the rank analysis and subspace 
108 for eigenvalue decomposition. If one or more eigenval- 
ues dominate the others, the channel subspace is estimated 
by computing the eigenvectors that correspond to the domi- 
nant eigenvalues. The eigenvectors spanning the channel 
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subspace are written to memory for further use in the 
channel subspace projection unit 109 where reduced rank 
channel parameter estimates are produced by computing the 
projection of the (N^-N^ N^) original channel estimates 
per estimation time interval onto the channel subspace, 
yielding N c; , reduced rank channel parameter estimates per 
estimation time interval. The results of the channel subspace 
projection unit 109 are written to memory for use in the 
demodulator 110. Optionally the channel subspace projec- 
tion unit 109 could generate equivalent full dimension 
channel parameter estimates, by re -transforming the N ch 
reduced rank channel parameter estimates into (N^-N^-N^) 
equivalent full dimension channel parameter estimates per 
estimation time interval. For example in a conventional 
RAKE-receiver design for the full rank channel model, the 
number of rake fingers for a full rank demodulator would be 
(Nj^-Nj^-Ne). A full rank demodulator would then use the 
(Nj^N^-Ne) original channel parameter estimates for the 
finger coefficients. Areduced complexity demodulator could 
eventually use only N ch RAKE fingers using the N ch 
reduced rank channel estimates as coefficients. However, 
since the receiver would generally be designed in anticipa- 
tion of a worst case situation, i.e., wherein (N^-N^-N^ 
fingers are implemented, it would be sufficient to compute 
(N^.-N^.-Nfr) correlated channel parameter estimates with 
improved estimation quality over of the N CA reduced rank 
channel parameter estimates. 

The rank analysis and subspace estimation unit 108 and 
the subspace projection unit 109 may be implemented in a 
Drgital Signal Processor (DSP), dedicated hardware, 
software, firmware, or a combination thereof. Modules 
wrthin receiver 100 may be incorporated together, and arc 
illustrated as separate blocks for clarity based on function. 

An exemplary configuration of one embodiment is illus- 
trated in FIG. 8 for a system having four (4) transmitter 
antennas and two (2) receiver antennas. Three (3) transmis- 
sion paths arc illustrated and labeled 1, 2 and 3. The points 
oi reflection for paths 1 and 2 are both on a same ellipse, 
wherein the ellipse is formed such that Tx and Rx are the 
focal points. Note that the ellipse is superimposed on the 
illustration of the physical layout of the system. Path 3 falls 
outsrde of the illustrated ellipse. Paths f and 2 have the same 
significant delay, t I , with respect to the receiver, while path 
3 has a significant delay, x 2 different from x 1 . The path delay 
rs a function of the configuration of the antennas as well as 
if the system. As illustrated, the four (4) 
is and the two (2) receiver antennas result 
in eight (8) channels. Each of the path delays, x 1 and %-,, 
pioduce an echo, wherein (N £ =2). The dimension oi the 
covariance matrix is given as (N^-N^..-N £ ) or sixteen (f 6) 
corresponding to the. (N^-N^-N^) channel rmpulse 
response samples. Therefore, the full rank channel param- 
eter vector is a 16-dimension vector. Using the rank reduc- 
tion methods described herein, the rank of the channel 
estimation may be reduced to three (3) dimensions, corre- 
sponding to paths 1, 2, and 3, wherein (N r ,=2). Note that 
where the mapping of the N c; , transmission paths to the 
(Nj-^.-N^.-Ne) channel impulse response samples is not 
known, the subspace may be extracted from configuration 
information. If the location and characteristics, such as 
direction and directionality, of the antennas are known, the 
information may be used to generate an arrav response or 
steering vector. Using the steering vector and path drrection 
information, which is also extractable using subspace 
algorithms, the angle of transmission, a, is estimated. If the 
antenna configuration has a fixed deployment the angle of 
calculable. A vector is formed including an 
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i for c 



20 



angle of tr; 

Similarly, an arrival angle vector is formed considering the 
receiver antennas. A linear transformation for the mapping 
of the N c „ transmission paths to the (N^N^-N^) channel 
impulse response samples is constructed using this informa- 5 
tion from both the transmitter and receiver configurations. 
This provides the matrix B as given in (19) hereinabove 
describing the linear transformation. The covariance matrix 
is derived therefrom as in (20) hereinabove. The process 
then proceeds as for the case where corresponding informa- 10 
tion is obtained from a priori knowledge. 

While one embodiment has been described herein with 
respect to the time domain, an alternate embodiment per- 
forms a rank reduction of the covariance matrix or a sample 
matrix in the irequency domain. If the parameters and 15 
equations are developed in the frequency domain, the pro- 
cess to estimate the channel then incorporates the frequency 
domain values. 

The previous description of the preferred embodiments is 
provided to enable any person skilled in the art to make or 20 
use the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the 
art, and the generic principles defined herein may be applied 
to other embodiments without the use of the inventive 
faculty. Thus, the present invention is not intended to be 25 
limited to the embodiments shown herein but is to be 
accorded the widest scope consistent with the principles and 
novel features disclosed herein. 

1. A wireless communication apparatus, comprising: 30 
a correlator operative to estimate a covariance matrix 

representing a link with a transmitter based on signals 
received from the transmitter; 

a rank analysis unit coupled to the correlator and opera- 
tive to estimate a rank of the covariance matrix; and " 

a channel estimation unit coupled to the rank analysis unit 
and operative to generate a reduced rank channel 

2. The apparatus of claim 1, wherein the covariance 4Q 
matrix represents a plurality of impulse responses between 
the apparatus and the transmitter. 

3. The apparatus of claim 1, wherein the correlator is 
operative to determine a correlation of at least two channels. 

4. The apparatus of claim 3, wherein the rank analysis unit 4J 
is operative to determine an eigenvalue corresponding to the 



7. A method in a wireless communication system, com- 
prising: 

estimating a covariance matrrx representing a link with a 
based on signals received from the trans- 



estrmatmg a rank of the covariance matrix; and 
generating a reduced rank channel estimate. 

8. An apparatus as in claim 1, wherein the correlator is 
operative to estimate delays in signals received. 

9. An apparatus as in claim 1, wherein the rank analysis 
unit is operative to compare the rank of the covariance 

10. An apparatus as in claim 9, wherein if the rank of the 
covariance matrix is less than the maximum rank, the 
channel estimation means is operative to generate a full rank 
channel estimate. 

11. A method as in claim 3, further comprising: 
determining if the rank of the covariance matrix is less 

than a maximum rank; and 
if the rank is less then NxM performii: 
a subspace of the matrix. 

12. A method as in claim 11, wherein e 
covariance matrix comprises: 

determining a matrix describing parametric relations of 
the link. 

13. A method as in claim 12, further comprising: 
deriving channel impulse responses for each channel 

based on the extracted subspace of the matrix; and 
demodulating a received signal using the channel impulse 
responses. 

14. A wireless communication apparatus, comprising: 



means for es 
link with < 

the i 

means for estimating a rank of the covariance i 
means for generating a reduced rank channel 
15. A wireless communication apparatus as ir 
further comprising: 

means for determining if the rank of the 
rank; and 



5. The apparatus of claim 4, wherein the rank analysis unit 
is operative to compare the estimated rank of the covariance 
matrix to a predetermined full value . s 

6. A wireless communication apparatus, comprising: 

a correlator means operative to estimate a covariance 

matrix representing a link with a transmitter based on 

signals received from the transmitter; 
a rank analysis means coupled to the correlator means and 5 

operative to estimate the rank of the covariance matrrx; 

and 

a channel estimation means coupled to the rank analysis 
means and operative to generate a reduced rank channel 



. means for performing an extraction of a subspace of the 
matrix if the rank is less then the maximum rank. 
16. A wireless communication apparatus as in claim 15, 
wherein the means for estimating a 
prises: 



17. A wireless communication apparatus as in claim 16, 
further comprising: 

means for deriving channel impulse responses for each 
j channel based on the extracted subspace of the matrix; 

means for demodulating a received signal using the 
channel rmpulse responses. 



